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ABSTRACT 

Polarization  adds  another  dimension  to  the  spatial  intensity  and  spectral  information  typically  acquired  in  remote 
sensing  Polarization  imparted  by  surface  reflections  contains  unique  and  discriminatory  signatures  which  may 
augment  spectral  target-detection  techniques.  While  efforts  have  been  made  toward  quantifying  the  polarimetric 
bidirectional  reflectance  distribution  function  (BRDF)  responsible  for  target  material  polarimetric  signatures, 
little  has  been  done  toward  developing  a  description  of  the  polarized  background  or  scene  clutter.  An  approach 
is  presented  for  measuring  the  polarimetric  BRDF  of  background  materials  such  as  vegetation. 

The  governing  equation  for  polarized  radiance  reaching  a  sensor  aperture  is  first  developed  and  serves  as  a 
basis  for  understanding  outdoor  polarimetric  BRDF  measurements,  as  well  as  polarimetric  remote  sensing.  The 
polarimetric  BRDF  measurements  are  acquired  through  an  imaging  technique  which  enables  derivation  of  the 
BRDF  variability  as  a  function  of  the  ground  separation  distance  (GSD).  An  image  subtraction  technique  is 
used  to  minimize  measurement  errors  resulting  from  the  partially  polarized  downwelled  sky  radiance.  Quan¬ 
tifying  the  GSD-dependent  BRDF  variability  is  critical  for  background  materials  which  are  typically  spatially 
inhomogeneous.  Preliminary  results  from  employing  the  measurement  technique  are  presented. 

Keywords:  Polarization,  BRDF,  clutter,  target  detection,  multispectral 

1.  INTRODUCTION 

Polarimetric  remote  sensing  in  the  visible  to  near-infrared  (VNIR),  or  in  the  spectral  regime  dominated  by  solar 
energy,  has  long  been  shown  to  offer  advantages  over  intensity-only  imaging.1-4  Benefits  such  as  improving  man¬ 
made  object  detection  are  often  touted,  as  well  as  possible  improvements  to  spectral  algorithms  used  for  detection 
and  identification.5  However,  virtually  all  efforts  fail  to  cast  polarimetric  remote  sensing  within  a  cohesive 
framework  in  which  a  priori  predictions  of  the  target  polarized  radiance  are  made,  as  is  done  with  spectral 
remote  sensing  techniques.  There  is  also  a  need  to  accurately  represent  and  model  the  polarized  background  or 
clutter  environment.  The  ability  to  model  target  and  background  polarimetric  signatures  is  a  prerequisite  for 
exploring  the  benefits  polarization  adds  to  spectral  target  detection  and  classification  algorithms.  Polarimetric 
signature  models  may  then  be  implemented  into  radiometrically  accurate  synthetic  image  generation  programs. 

First,  the  bidirectional  reflectance  distribution  function  (BRDF)  is  introduced  and  the  more  general  po¬ 
larimetric  BRDF,  which  quantifies  polarized  signatures.  Next,  the  governing  equation  for  polarized  radiance 
reaching  a  sensor  is  developed.  The  equation  highlights  the  role  of  the  polarimetric  BRDF,  and  also  serves  as  a 
foundation  for  polarimetric  remote  sensing.  The  governing  equation  guides  the  implementation  of  the  outdoor 
measurement  technique  for  quantifying  background  material  polarimetric  signatures.  The  imaging  approach  for 
these  measurements  is  reviewed,  and  characterizations  of  the  system  described.  Finally  preliminary  results  are 
shown. 
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2.  BIDIRECTIONAL  REFLECTANCE  DISTRIBUTION  FUNCTION 

The  bidirectional  reflectance  distribution  function  (BRDF)  provides  a  complete  description  of  the  reflectance 
properties  of  a  material.  Specifically,  BRDF  quantifies  the  geometric  distribution  of  radiance  reflected  from  a 
surface  illuminated  by  a  source  in  an  arbitrary  position  above  the  hemisphere  of  the  material. 


The  BRDF  is  given  by 


dLr(Or,  (for') 
dE(0U<f>i) 


(1) 


where  Lr  is  the  surface  leaving  spectral  radiance  [55^5]  and  E  is  the  spectral  irradiance  [^] ,  resulting  in 

BRDF  having  units  of  sr-1.*  The  incident  and  reflecting  zenith  angles  are  8i,0r  and  the  corresponding  azimuth 
angles  <ph<f>r .  The  nomenclature  used  here  is  that  recommended  by  Nicodemus,6’ 7  which  has  subsequently  been 
adopted  by  many  authors. 


The  most  general  expression  for  BRDF  fully  incorporates  all  incident  and  reflected  polarization  states  as 
well.8,9  This  expression  may  be  cast  as 


Lr(6r ,  cj)r )  —  Fr(0j,  4>i ;  0r,  4>r\ 
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where  now  the  BRDF  (Fr)  is  a  4  x  4  Mueller  matrix  and  the  incident  irradiance  and  surface  leaving  radiance 
are  given  as  Stokes  vectors,  E  and  Lr. 

In  passive,  VNIR  remote  sensing,  the  polarization  ellipticity  may  be  considered  negligible,1’10  and  therefore 
the  4  x  4  BRDF  Mueller  matrix  may  be  reduced  to  a  3  x  3  representation.  Therefore,  (2)  is  given  explicitly  by 


Lo  ‘ 

foo  foi  fo2 

Eq 

Li 

= 

flO  fll  fl2 

Ex 

l2 

1 - 

0 

S- 

to 

E2 

and  the  foo  component  of  Fr  is  seen  to  equal  the  scalar  BRDF  value,  fr  of  (1). 

Many  man-made  materials,  which  often  comprise  “targets”  in  spectral  detection  algorithms,  typically  have 
spatially  uniform  surfaces  over  which  the  BRDF  does  not  vary  significantly  (e.g.,  painted  surfaces).  This  results 
in  resolved  target  surfaces  of  similar  orientation  having  minimal  pixel-to-pixel  variance.  Laboratory  BRDF 
measurements  successfully  characterize  such  surfaces11’12  and  these  empirical  measurements  are  used  to  fit 
BRDF  models,13-15  enabling  a  priori  predictions  of  spectral  radiance  received  by  a  sensor  given  arbitrary  solar 
and  sensor  orientations. 

Similarly,  materials  comprising  the  “background,”  which  are  usually  vegetation  and  soils,  have  been  studied 
extensively  for  their  anisotropic  reflectance  characteristics.15  The  BRDF  characterization  of  these  surfaces  are 
used  toward  deducting  biophysical  parameters  for  agricultural  applications  or  surface  albedo  calculations  for 
global  climate  considerations.17  The  spatial  extent  over  which  these  measurements  are  made  is  usually  orders 
of  magnitude  greater  than  that  of  target  material  laboratory  measurements.  This  is  motivated  in  part  by  the 
relatively  large  Ground  Separation  Distance  (GSD)  of  earth  resource  satellites  which  monitor  agricultural  and 
global  climate  processes. 

Characterizing  the  BRDF  of  background  materials  at  a  comparable  spatial  extent  to  that  of  target  materials 
is  hampered  by  the  spatial  variability,  or  texture  of  natural  materials.  The  background  has  significant  signature 
variability  at  GSDs  of  interest  for  target  detection  applications.  This  within-material  BRDF  variability  has  been 
termed  the  bidirectional  texture  function18  or  the  bidirectional  reflectance  variance  function  (BRVF).19  BRVF 
may  be  considered  a  BRDF  probability  distribution  function  for  a  given  GSD.  As  the  GSD  increases,  the  BRDF 
variability  decreases  since  more  of  the  texture  is  averaged  within  a  single  pixel. 

The  variability  exhibited  by  background  materials  motivates  the  need  for  a  measurement  approach  which  cap¬ 
tures  the  variability.  An  imaging  technique  has  been  previously  reported  by  the  authors  as  a  suitable  method,20 
and  it  is  the  polarimetric  implementation  of  this  approach  which  will  be  discussed. 


'All  radiometric  quantities  are  assumed  to  be  spectral,  though  not  always  explicitly  shown. 


3.  GOVERNING  POLARIZED  RADIANCE  EQUATION 

Prior  to  discussing  the  polarimetric  BRDF  measurement  approach,  it  is  instructive  to  first  derive  an  expression 
for  the  polarized  radiance  reaching  a  sensor  aperture.  This  establishes  the  framework  from  which  quantitative 
polarimetric  remote  sensing  must  operate,  and  also  guides  the  measurement  technique.  The  radiometric  equation 
for  the  unpolarized  intensity  is  first  introduced  and  nomenclature  established.  The  polarimetric  representation  of 
these  equations  are  then  derived,  which  results  in  emphasizing  the  role  of  the  polarimetric  BRDF.  The  polarized 
governing  equation  is  not  known  to  have  been  previously  addressed  in  detail. 

3.1.  Governing  Equation — Magnitude 

The  total  radiance  in  the  visible  to  near  infrared  (VNIR)  portion  of  the  spectrum  (i.e.  that  of  solar  origin) 
reaching  a  sensor  aperture  ( Ls )  may  be  approximated  as  the  sum  of  three  radiance  sources: 

1.  direct  solar  reflections  from  the  target,  Lr 

2.  target-reflected  downwelled  radiance  from  the  skydome,  Ld 

3.  upwelled  atmosphere  radiance  resulting  from  solar  scatter  in  the  atmosphere  along  the  target  to  sensor 
path ,  La 

The  order  of  the  radiance  terms  above  is  that  of  typically  decreasing  magnitude,  though  the  ground  or  target 
reflectance  and  atmospheric  conditions  greatly  influence  their  relative  values  (c.f.  Fig.  4.12,  Tbl.  4.1  of21).  These 
radiance  terms  are  functions  of  the  incident  and  reflected  zenith  and  azimuth  angles,  0'T  and  azimuth  angle 

4>  ■ 

An  expression  for  the  radiance  from  the  direct  solar  reflection,  Lr,  is  obtained  by  first  considering  the 
exoatmospheric  solar  irradiance,  Es,  which  propagates  through  the  atmosphere  along  the  solar-to-target  path 
having  a  transmittance  of  Ti .  When  incident  upon  a  surface,  it  is  then  reflected,  and  again  attenuated  by  the 
atmosphere  along  the  ground-to-sensor  atmospheric  path  by  rr.  Often  the  reflectance  is  considered  Lambertian, 
or  isotropic,  such  that  a  reflectance  factor  is  used  as  an  approximation  to  what  is  properly  the  bidirectional 
reflectance  distribution  function  (BRDF),  fr.  Assembling  these  terms,  Lr  may  therefore  be  expressed  as 

Lr  =  Tr(0'r)  fr(Oi,0r,4> )  COS#)  Ti(#))  £s(#))  ,  (4) 

where  the  unprimed  coordinate  system  is  one  relative  to  the  normal  of  a  material  surface. 

In  a  similar  fashion,  target-reflected  radiance  from  the  sky,  Ld  may  be  derived.  The  downwelled  radiance 
distributed  over  entire  hemisphere  of  the  sky,  L^\  is  integrated  to  sum  irradiance  contributions  onto  the  target 
across  the  sky,  which  is  modified  by  the  cosine  of  the  incident  angle  upon  the  surface  normal.  As  before, 
each  of  these  irradiance  contributions  is  then  reflected  by  the  surface  BRDF,  which  is  then  attenuated  by  the 
target-to-sensor  atmospheric  transmittance  as  before.  Replacing  the  BRDF  by  an  isotropic  reflectance  factor 
greatly  simplifies  the  expression,  as  the  reflectance  factor  may  be  placed  outside  the  integral.  However,  the  more 
stringent  BRDF  must  be  retained  as  it  is  essential  to  polarimetry.  An  appropriate  expression  for  Ld  is  therefore 

Ld  =  Tr(0'r)  J f  MOi,  0r,  4>)  cos  #'  L%< (#),  t')  dnt ,  (5) 

fit 


where  dfi)  =  sin  d6{d4>  . 

A  representation  for  the  upwelled  atmospheric  radiance,  Lu  will  not  be  attempted,  as  it  is  rather  complex  and 
usually  approximated  by  atmospheric  scattering  codes  such  as  MODTRAN,22  as  is  the  downwelled  sky  radiance 
component  (L% ‘ )  in  equation  5.  The  upwelled  radiance  is  given  simply  in  order  to  show  the  geometry  dependence 
as 

Lu  =  Lu(er,cf>)  .  (6) 


3.2.  Governing  Equation — Stokes  Representation 

Transforming  (4-6)  into  the  polarized  representation  is  accomplished  using  the  Mueller-Stokes  formalism  com¬ 
monly  used  in  polarized  radiometry.  In  brief,  all  radiometric  flux  values  are  replaced  by  Stokes  vectors  and 
“transfer”  functions  such  as  atmospheric  transmittance  and  reflectance  (BRDF)  are  replaced  by  Mueller  matri- 
ces. 

Prior  to  making  these  substitutions,  some  simplifications  are  appropriate.  First,  the  exoatmospheric  solar 
irradiance  may  be  considered  randomly  polarized,  so  only  the  scalar  magnitude  (or  first  Stokes  component)  of  the 
direct  solar  irradiance  need  be  considered.  Second,  the  atmospheric  transmittance  values  in  (4-6)  all  represent 
forward  scattering,  which  retains  the  incident  polarization.  Therefore,  the  scalar  values  for  t*  and  rr  may  be 
used  without  having  to  resort  to  a  Mueller  matrix  representation.  Equations  4-6  therefore  become 
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where  Fr  is  now  the  polarimetric  BRDF  (pBRDF).  Some  knowledge  of  the  upwelled  polarized  radiance  ( Lu ) 
along  the  target  and  sensor  may  be  gained  from  Rayleigh  scattering  theory  and  other  sources  such  as  Coulson 
and  Chandrasekhar.24  However,  knowledge  of  the  polarized  downwelled  radiance,  Ld  1  is  more  problematic  since 
this  term  often  has  a  high  spatial  variability,  e.g .,  varying  cloud  cover. 

The  total  polarized  radiance  reaching  a  sensor  aperture  is  then 

Ls  —  Lr  +  Ld  +  Lu  .  (1^) 

Atmospheric  scattering,  generally  proportional  to  A-4,  results  in  Ld  and  Lu  having  relatively  large  magnitudes  at 
shorter  wavelengths  compared  to  Lr)  especially  from  orbital  altitudes.  Atmospheric  polarimetric  remote  sensing 
uses  this  phenomena  to  minimize  ground  reflected  polarization  signatures  to  better  extract  atmospheric  water 
vapor  and  aerosol  properties.25 

Similarly,  for  polarimetric  remote  sensing  of  land  features  one  wants  the  magnitude  of  the  direct  solar  reflected 
radiance  to  be  large  compared  to  the  reflected  radiance  from  the  downwelled  sky  and  upwelled  atmospheric 
scattering,  i.e Lr  >  Ld,Lu.  This  provides  optimal  conditions  for  estimating  the  polarimetric  BRDF,  Fr. 
Exploiting  polarimetric  signatures  in  a  manner  analogous  to  spectral  signatures  requires  estimating  Fr  given  the 
polarized  radiance  reaching  the  aperture,  Ls. 

Estimating  Fr  given  the  radiance  at  the  sensor  aperture  proceeds  as 

Lr  —  Ls  —  Ld  Lu  (11) 

rr  Fr  Ti  cos  9i  Es  =  Ls  —  rr  j  j  Fr  cos 6^  Ld%  dfti  —  Lu  •  (12) 

ft* 

Since  the  exoatmospheric  irradiance  is  randomly  polarized,  only  the  first  column  of  the  pBRDF  Mueller  matrix 
is  of  concern  in  the  Lr  expression.  In  fact,  overhead  polarimetric  remote  sensing  may  only  retrieve  the  first 
column  of  the  polarimetric  BRDF  matrix.  (Solving  for  other  matrix  elements  requires  illumination  by  varying 
polarization  states).  With  this  consideration  (12)  may  be  expressed  as 

Ti  cos  O'i  Es  =  Ls  -  rr  JJ  Fr  cos^  Ldl  -  Lu  (13) 

ft* 

with  the  first  column  of  the  pBRDF  given  by 

Ls-TrJ fQi  Fr  COS  dtij  -  L, 
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Solving  for  Fr  is  complicated  by  its  inclusion  in  the  integral  of  the  Ld  term,  which  also  contains  the  highly 
spatially  variable  and  generally  unknown  downwelled  radiance  component,  Ldf.  However,  under  nominal  sky 
conditions,  the  magnitude  of  the  direct  solar  irradiance  for  A  >  600  nm  is  five  times  that  of  the  integrated 
sky-dome  irradiance,  increasing  to  lOx  for  A  >  1000  nm.  This  makes  it  reasonable  to  approximate  the  polarized 
radiance  contribution  of  the  downwelled  sky  radiance  as  an  error  term. 
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Therefore,  polarimetric  remote  sensing  may  recover  the  first  column  of  polarimetric  BRDF  Mueller  matrix  to 
within  the  error  resulting  from  downwelled  sky  radiance,  presented  as 


/oo  +  e0 
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/20  +  e2 
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Note  that  €0  is  always  positive,  and  for  diffuse  surfaces  the  ratio  of  is  equivalent  to  the  ratio  of  the  downwelled 
sky  irradiance  to  the  direct  solar  irradiance.  The  linear  polarization  terms  ex  and  e2  may  either  be  positive  or 
negative  and  represent  the  influence  of  the  partially  polarized  downwelled  sky  radiance. 


4.  POLARIMETRIC  BRDF  MEASUREMENT  APPROACH 

Ideally,  BRDF  measurements  are  made  in  a  lab  environment  using  a  “point”  illumination  source  with  careful 
control  and  minimization  of  stray  light.  However,  many  materials  such  as  vegetation  do  not  lend  themselves  to 
easy  indoor  measurements  due  to  alteration  of  their  natural  state  or  simply  from  their  physical  size  ( e.g .,  a  tree 
canopy).  Outdoor  BRDF  measurements  of  such  materials  becomes  a  necessity,  and  many  approaches  have  been 
successfully  employed.26-29  Wide  field  of  view  (FOV)  imaging  systems  may  be  used  which  efficiently  enable  the 
simultaneous  measurement  of  multiple  scattering  angles.30-32 

Our  approach  is  having  a  narrow  FOV  («  10°)  imaging  system  to  make  BRDF  measurements.  Each  image 
pixel  is  approximated  as  the  same  scattering  angle  as  that  at  the  center  of  the  image,  such  that  the  average 
radiance  across  the  focal  plane  enables  determination  of  the  BRDF.  Such  an  approach  limits  the  scattering  angle 
resolution  to  the  FOV,  but  this  is  not  a  concern  for  most  natural  surfaces  which  are  not  appreciably  specular 
and  hence  do  not  have  rapid  BRDF  changes  over  the  10°  FOV  of  the  system.  The  impetus  for  this  technique  is 
the  ability  to  quantify  the  the  BRDF  variability  (as  discussed  in  §2).  Multiple  scattering  angles  are  sampled  by 
repositioning  the  camera  in  the  hemisphere  above  the  measurement  surface. 

This  technique  may  be  used  at  any  distance  from  the  measurement  surface — the  only  prerequisite  is  that 
the  ground  FOV  (GFOV)  is  large  enough  that  it  adequately  integrates  the  spatial  variability  or  texture  of  the 
material.  For  instance,  a  GFOV  of  1  foot  may  be  adequate  for  grass,  asphalt  and  aggregate,  but  measurements 
of  tree  canopies  and  shrubs  would  require  a  larger  GFOV.  For  easy  field  use  not  requiring  elevated  platforms 
or  other  positioning  devices,  an  operating  distance  for  the  measurements  discussed  here  was  6  feet,  providing  a 
GFOV  of  approximately  1  foot. 


4.1.  BRDF  Measurement 

A  successful  technique  for  outdoor  BRDF  measurements  may  be  developed  by  considering  the  radiance  contri¬ 
butions  to  a  sensor  (c./.  equation  (10)).  It  is  first  noted  that  imaging  surfaces  at  a  distance  of  6  feet  results  in 
negligible  atmospheric  scattering  along  the  surface  to  sensor  path,  such  that  Lu  ~_0.  The  surface  radiance  is 
therefore  composed  of  the  direct  solar  and  downwelled  sky  reflectance,  or  Lr  and  Ld-  The  measurement  made 
when  the  surface  is  illuminated  by  the  sun  and  downwelled  sky  radiance  will  be  referred  to  as  image  C  (Figure 
2). 


The  downwelled  sky  radiance  is  a  stray  light  source  for  the  purpose  of  BRDF  measurements.  It  may  be  directly 
measured  and  eliminated  via  an  image  subtraction  technique.  La  is  measured  by  placing  the  measurement  area 
in  shadow,  and  imaging  the  shadowed  surface  (image  D).  In  this  manner  it  is  seen  that 

LrocC-D  =  (Lr  +  Ld)-(Ld)  .  (18) 


The  error  terms  shown  in  equation  (17)  are  therefore  eliminated  by  the  “shadow”  image.  This  is  quite 
valuable,  as  comparison  of  the  C  and  C  —  D  data  quantifies  the  change  to  the  linear  Stokes  components  resulting 
from  the  sky  polarization. 


4.2.  Radiance  Calibration 

The  “digital  counts”  recorded  by  the  imaging  system  may  be  transformed  into  absolute  radiance  levels  by  use 
of  a  Spectralon  calibration  target.  Spectralon  has  a  highly  Lambertian,  angular-invariant  BRDF  of  £,  with  a 
randomly  polarized  reflectance  of  p  >  0.97  across  most  of  the  VNIR  spectrum.33  As  with  the  surface  being 
measured,  images  of  the  calibration  target  are  taken  both  in  sun  and  in  shadow,  images  A  and  B,  respectively. 

When  acquiring  multiple  images  over  a  short  time  period  such  that  the  atmospheric  conditions  and  solar  zenith 
position  (0i)  do  not  change  appreciably,  the  BRDF  may  be  determined  by  the  ratio  of  the  known  calibration 
target  BRDF  to  that  of  the  unknown  surface  or 
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In  terms  of  the  digital  counts  of  the  pixels  in  each  of  the  four  images,  A  through  D,  the  BRDF  is 

pTC-D' 

Jr  7T  A-B 


(20) 


When  imaging  a  calibration  target  such  that  it  occupies  the  full  FOV,  this  technique  also  self-corrects  for  the  so 
called  “lens  falloff”  irradiance  reduction  away  from  the  center  of  the  focal  plane. 


4.3.  Polarimetric  BRDF 

The  polarized  radiance  leaving  the  surface  may  be  quantified  as  a  Stokes  vector  using  well-established  ap¬ 
proaches.8,34  In  our  implementation,  images  of  the  surface  are  acquired  under  four  different  linear  polarization 
filter  orientations  relative  to  the  horizon:  0°,  45°,  90°  and  135°.  This  enables  derivation  of  the  Stokes  vector 
according  to 

|  Sn  I  I  | 

(21) 

where  J**  represents  an  image  acquired  with  the  polarization  filter  set  at  xx°.  It  is  noted  that  the  first  Stokes 
component  is  derived  using  an  average  of  both  sets  of  cross-polarized  images. 

In  terms  of  the  images  using  the  calibration  target,  it  is  seen  from  equations  (20,  21)  that  the  polarimetric 
BRDF  is  therefore 
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where  arb  indicates  an  arbitrary  polarization  filter  orientation  for  imaging  the  calibration  target,  since  this 
radiance  is  randomly  polarized. 

To  summarize,  for  each  hemispherical  scattering  position,  a  total  of  8  images  are  acquired  of  the  target  surface, 
4  polarization  orientations  with  2  illumination  conditions  (full  sun  and  shadow).  A  minimum  of  2  calibration 
target  images  must  be  taken,  one  for  each  illumination  condition.  Therefore  a  data  set  at  one  scattering  position 
comprises  10  images. 


RGB  Histograms 


Digital  Count 


Figure  1.  The  RGB  BRDF  distributions  or  BRFV  for  a  “grass”  measurement.  Histograms  are  shown  for  the  full  image 
resolution  and  at  an  arbitrarily  defined  “GSD  =  1”.  The  averaging  of  the  texture  as  a  function  of  GSD  is  illustrated  at 
right. 


4.4.  BRDF  Probability  Distribution  (BRVF)  Calculation 

Thus  far,  only  the  average  digital  count  values  over  the  entire  image  have  been  considered  in  deriving  the  BRDF. 
However,  the  impetus  for  this  technique  is  the  ability  to  quantify  the  BRDF  variability,  or  BRVF  discussed  in 
§2.  The  variability  is  obviously  a  function  of  the  GSD,  as  a  larger  GSD  results  in  greater  averaging  of  texture 
within  a  pixel,  and  hence  decreased  pixel  to  pixel  variability. 

The  high-resolution  images  acquired  with  the  BRDF  measurement  system  may  then  be  used  to  generate 
the  BRVF  given  the  anticipated  GSD  of  a  remote  sensing  sensor.  Generating  the  BRVF  is  accomplished  by 
convolving  the  image,  f[x,y],  with  a  convolution  kernel  h[x,y]  representative  of  the  GSD  of  interest.  The  result 
is  a  low-pass  filtered  image,  g[x,  y\  with  the  spatial  texture  representative  of  the  GSD  of  h[x,  j/].  This  is  presented 
mathematically  as 

^  771  —  1  71  —  1 

K  y\  =  h^x  -  *> y  -  ft  (23) 

^  i- 0  j— 0 

where  X 2  is  a  weighting  factor  such  that  the  average  magnitude  of  g[x,y]  is  that  of  the  original  image,  f[x,y]. 
Ideally  h[x,y ]  is  the  point  spread  function  of  the  remote  sensing  platform  in  question,  but  for  quick  processing 
a  simple  function  with  a  unit  magnitude  and  square  spatial  extent  is  used  (termed  a  RECT  function  by  some35). 
Figure  1  illustrates  the  effect  using  a  simple  color  (RGB)  image  of  “grass”  taken  with  a  commercial  digital 
camera. 

Unlike  the  polarimetric  BRDF  determination,  the  accuracy  of  the  BRVF  depends  upon  the  degree  of  the 
spatial  registration  of  the  four  sets  of  polarized  C  and  D  images.  When  the  size  of  the  convolution  kernel  is 
commensurate  with  the  spatial  registration  accuracy,  significant  errors  result.  The  same  is  true  of  movement 
of  measurement  surface  while  acquiring  the  four  polarization  orientations,  e.g.  grass  blowing  in  the  wind.  As 
required,  the  C  and  D  image  sets  should  be  spatially  registered  prior  to  performing  BRVF  calculations. 

A  summary  of  the  general  measurement  steps  for  this  technique  is  presented  as  Figure  2.  Depending  on  the 
polarimetric  imaging  system  used  to  make  the  measurements,  this  process  should  be  modified  accordingly,  such 
as  spectral  filter  changes,  etc. 


1.  Fix  Exposure  /  Gain  Settings 


2.  Radiance  Calibration 


3.  Surface  Measurement 


Figure  2.  Measurement  and  process  flow  for  making  polarimetric  BRDF  and  BRVF  measurements  with  the  camera 
system. 


Figure  3.  The  assembled  imaging  system  mounted  on  a  tripod. 


5.  IMAGING  SYSTEM  DESCRIPTION  &  CHARACTERIZATION 

The  imaging  system  consists  of  a  SenSys  1602E  camera  having  a  1536  x  1024  thermoelectric-cooled,  12-bit  silicon 
CCD  with  a  response  nonlinearity  <  0.5%.  A  filter  wheel  located  between  the  lens  and  the  CCD  is  used  to  mount 
25  mm  diameter  band-pass  filters.  The  spectral  filter  wheel  housing  accepts  a  standard  F-mount  lens,  to  which  a 
Nikon  50  mm,  f/1.8  lens  is  used.  A  linear  polarization  filter  is  mounted  external  to  the  lens  on  an  optics  post  in 
a  precision  rotary  mount,  which  is  mounted  to  a  common  optics  board  with  the  camera.  This  assembly  is  then 
mounted  on  a  tripod.  To  demonstrate  the  technique,  data  is  only  taken  at  two  spectral  bands,  550  ±  5  nm  and 
750  ±  12  nm.  An  overview  of  the  system  is  shown  in  Figure  3. 

The  imaging  system  was  characterized  in  order  to  gain  an  understanding  of  the  measurement  uncertainties 


and  limitations.  First,  it  was  noted  the  “dark”  images  of  the  camera  were  highly  repeatable,  and  have  negligible 
error  contribution  to  the  series  of  images  used  to  make  a  measurements. 

The  lens  falloff,  or  focal  plane  irradiance  decrease  away  from  the  center  of  the  array,  was  also  quantified  by 
imaging  into  an  integrating  sphere,  which  provided  a  uniform  radiance  field.  At  an  aperture  setting  of  f/8.0, 
where  the  system  is  usually  operated,  the  edge  of  the  measurement  area  on  the  focal  plane  has  a  response  of 
0.94  ±  0.01  relative  to  that  at  the  center  of  the  array.  Correction  to  the  lens  falloff  is  only  necessary  under 
circumstances  when  the  calibration  target  may  not  be  imaged  over  the  full  FOV  of  the  system. 

5.1.  Polarization  Filter  Alignment  Errors 

A  thorough  error  analysis  of  the  polarization  errors  is  quite  involved  and  beyond  the  scope  of  this  paper.  Only 
the  final  error  equations,  in  terms  of  the  Stokes  components,  and  the  anticipated  net  error  are  presented.  First, 
there  is  uncertainty  in  the  polarization  filter  orientation  relative  to  the  local  horizon,  which  is  considered  the 
absolute  reference  frame  for  the  polarization  orientation.  This  error  is  not  as  critical,  as  it  does  not  impact  the 
measured  DOP,  but  only  the  relative  magnitude  between  the  Si  and  52  Stokes  components. 

More  important  are  the  errors  resulting  from  manually  positioning  the  polarization  filter  in  the  four  orien¬ 
tations.  This  error  is  complex  and  is  a  function  of  the  incident  polarization  magnitude  and  direction,  and  of 
course  the  polarization  filter  orientation  error,  exx  for  the  xx°  filter  position.  It  may  be  shown  that  the  error  in 
the  derived  intensity  or  first  Stokes  component  is  given  by 

So  =  50in  4-  e0tot  =  50in  +  ^p{[cos(2e0)  -  cos(2e9o)]  +  [sin(2ei35)  -  sin(2e45)]} 

+  ^L|[cos(2e45)  -  cos(2e135)]  +  [sin(2e0)  -  sin(2e9o)]| 

with  errors  in  the  linear  Stokes  components  given  by 

51  =  SUn  +  eitot  =  i|5iin  [cos(2e0)  +  cos(2e90)]  +  52in  [sin(2e0)  +  sin(2e90)]| 

52  =  S2in  ±  e2tot  =  ^52in  [cos(2e45)  -I-  cos(2ei35)]  -  Slin  [sin(2e45)  +  sin(2ei35)]} 

It  is  estimated  that  the  positioning  accuracy  of  the  polarization  filters  is  exx  —  0  ±  0.25°  at  a 
level.  Numerical  simulations  of  equations  (24)-(26)  result  in  an  anticipated  2a  error  less  than 
Stokes  component. 

5.2.  Finite  Filter  Contrast  Error 

Next,  the  error  resulting  from  the  finite  contrast  or  cross  polarized  “leakage”  are  presented.  This  performance 
metric  may  be  given  by  the  cross-polarized  transmittance  relative  to  the  like-polarized  transmittance,  t®.  This 
error  produces  the  intuitive  result  of  overestimating  the  total  intensity,  50,  and  underestimating  the  linear  Stokes 
components.  Again  without  derivation,  this  results  in 

So  —  (1  +  T(g))So<n ,  Si  =  (1  —  T®)Si<n  and  =  (1  —  inj  (^) 

with  the  impact  on  DOP  given  by 

DOP  =  ( DOPin  .  (28) 

\1  +  W 

It  is  noted  that  this  error  is  systematic  and  may  be  corrected  with  knowledge  of  t®.  For  our  system,  r®  at  550 
nm  is  estimated  as  0.015. 


(24) 

(25) 

(26) 

2 a  confidence 
±1.2%  for  all 


Figure  4.  The  Stokes  and  DOP  images  of  a  “Magic  8-ball”  under  ambient  lighting  conditions  at  550  nm.  So  (top  left), 
Si  ( top  right),  S2  ( bottom  left)  and  DOP  ( bottom  right). 


5.3.  Test  Image 

Having  gained  an  understanding  of  the  system  performance,  a  data  set  was  acquired  by  imaging  a  “Magic  8-ball  ” 
The  ball  is  well-suited  for  demonstrating  polarization  phenomenology  as  it  has  a  highly  smooth,  specular  surface, 
including  regions  of  black  and  white  having  very  low  and  very  high  diffuse  (randomly  polarized)  reflectance.  In 
addition,  the  curvature  of  the  ball  provides  multiple  specular  view  angles.  The  ball  was  imaged  under  ambient 
lighting  conditions  in  front  of  a  Spectralon  panel.  The  images  were  processed  according  to  equation  (21)  providing 
the  Stokes  vectors,  from  which  the  DOP  was  calculated  (Figure  4).  The  DOP  image  provides  a  pleasing  result- 
reflectance  from  the  Spectralon  panel  off  the  edges  of  the  ball  provide  a  DOP  commensurate  with  that  expected 
from  Fresnel  reflectance,  with  a  peak  magnitude  reached  near  Brewster’s  angle. 

6.  PRELIMINARY  RESULTS 

J  plan  on  including  some  lab  measurements  of  “gravel”  to  demonstrate  the  anticipated  outdoor  results.  These 
indoor  measurements  will  admittedly  be  approximations }  due  to  the  lack  of  a  uniform  irradiance  source.  Will 
include  a  figure  similar  to  the  “8-ball”  one ,  and  also  generate  BRVF  statistics. 

7.  CONCLUSIONS 

The  polarimetric  governing  equation  for  reflected  radiance  has  been  developed  and  shown  to  provide  a  quan¬ 
titative  framework  from  which  polarimetric  remote  sensing  must  operate.  Measuring  background  material  po¬ 
larization  signatures  quantifies  the  noise  floor  for  polarimetric  target  detection  and  identification  algorithms. 
The  presented  measurement  technique  enables  polarized  signatures  tailored  to  the  point  spread  function  or 
GSD  of  specific  imaging  systems.  Such  measurements  may  be  applied  to  polarimetric  BRDF  models,  enabling 
radiometrically-accurate  synthetic  image  generation.  Synthetic  hyperspectral  polarimetric  imagery  allows  exten¬ 
sive  modelling  of  the  varying  conditions  under  which  spectral  and/or  polarimetric  target  detection  and  identifi¬ 
cation  algorithms  may  operate. 
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